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We review current experimental results on charm mixing and CP violation. We survey 
experimental techniques, including time-dependent, time-independent, and quantum- 
correlated measurements. We review techniques that use a slow pion tag from D*+ — > 
tt+D'' -|- c.c. decays and those that do not, and cover two-body and multi-body de- 
cay modes. We provide a summary of D-mixing results to date and comment on future 
experimental prospects at the LHC and other new or planned facilities. 
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1. Introduction 

Quantum-mechanical mixing between neutral meson particle and anti-particle flavor 
eigenstates provides important information about electroweak interactions and the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix, as well as the virtual particles that 
are exchanged in the mixing process itself. The two parameters characterizing D'^-D^ 
mixing are 

x=^^, AM = Ml- M2 (1) 

Ar 

y^^, Ar^ri-r2 (2) 
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where Mi. 2 are the masses of -D1.2, ri_2 are the decay widths, and F = (Fi + F2)/2 
is the mean decay width. 

Mixing between the states and ^ and 5°, and Bg and Bs is well es- 
tablished. Mixing in these systems is well described by standard model (SM) box 
diagrams containing up-type (u, c, t) quarks. In contrast, the D^-D'^ SM mixing am- 
plitude at short distances involves loops containing down-type (d, s, h) quarks. The 
s and d box amplitudeJl^ together are suppressed by by a factor (m^— m^)^/ (rri^ni^) 
due to the Glashow-Iliopoulos-Maiani (GIM) mechanism, while the contribution 
from loops involving b quarks is further suppressed by the Cabibbo-Kobayashi- 
Maskawa (CKM) factor iKbKtl ViKsKsP = 0(10"^). The contribution of the 
box diagrams shown in Fig.[Ha-b) to x is a;*""^ sa 5 x 10~^ [ms/0.2 GeV/c^] ^.121 The 

di-penguin diagram shown in Fig.[IJc) contributes at a similar level, but with oppo- 

9 . . 

site sign. . Such diagrams contribute only to x. A perturbative QCD next-to-leading 

order (NLO) analysis of D°-D° mixing using an operator product expansior l'^ l ^ ^ ^ l to 

evaluate AF in terms of local | AC| = 1 operators, followed by a dispersion relation 

to evaluate AM,[3 obtains: x,y ~6 x 10~^. Taken together, the short-distance SM 

predictions are x ~ 0(10"^), y ~ 0(10^''), far below the current measurements, 

x,y^OilO-'). 

The long-distance contributions to D'^-D^ mixing are inherently nonperturbative 
and thus difficult to estimate. There are two approaches to estimating D'^-D^ mixing 
in the SM: An inclusive approach uses the Operator Product Ex pansio n and quark- 
hadron duality to expand x and y in terms of local operatorsl^EEl jf ^j^g charm 
quark mass rric is large compared to the scale of strong interactions, the series can 
be truncated after a few terms. Such calculations typically yield x,y < 10^^. The 
exclusive approach sums over intermediate hadronic states to which both and 
D'^ can decay, as shown schematically in Fig. [2l Here, non- vanishing y arises from 
SU{3)f breaking in decay rates when summing over intermediate states within 
an SU{3)f multiplet. Ref. |8] found that SU{3)f violation in the final state phase 
space could provide enough SU{5)f breaking to generate y ^ 10~^. Ref. |7|used a 
dispersion relation to relate x to y and found x to be in the range (—1.0, —0.1) x y. 

New physics (NP) processes, some examples of which are shown in Figs, [ijd-f), 
could enhance the D^-D'^ mixing rate to the level of experimenta l detecti on, but the 
predictions for these rates also span many orders of mag nitude.l^li^. Given the 
uncertainties in both the SM and NP calculations, observation of D^-D^ mixing at 
0(10"^) does not unambiguously indicate the presence of new physics. See Refs.|9| 
and 1101 for a summary of D'^-D^ mixing parameter predictions. 

Evidence for D^-D^ mi xing w as reported in 2007 using high-luminosity data sets 
acquired at the B factorieJI^CSl and Tevatron collider. W While the significance of 
the current wo rld average for D'^-D'^ mixing is greater than ten standard devia- 
tions (10cr),flS! date no one single D'^-D'^ mixing measurement exceeds 5cr, the 
commonly accepted criterion for observation. 
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Fig. 1. Possible short-distance amplitudes contributing to D^-Z?" mixing, (a— b) SM boxes; (c) 
SM di-penguin; (d): new physics flavor-changing neutral current process mediated by a heavy Z'"; 
(e-f): charged Higgs in the mixing loop. 




Fig. 2. Long-distance contribution from an intermediate state / to D^-D^ mixing. H is the 
Hamiltonian governing weak decays. From Ref. 1161 



1.1. D°-D° Mixing Formalism 

The D'^ and 13° mesons are produced as flavor eigenstate with charm quantum 
numbers C — +1 and —1, respectively. They propagate and decay according to the 
Schrodinger equation: 



Mixing between D° and D° occurs because these flavor states are not the eigenstates 
Di and D2 of the D^-D'^ mass matrix M — ?r/2, but linear combinations of them. 
Assuming that the product of charge conjugation, parity a nd tirn e reversal (CPT) 
is conserved,!!^ the eigenstates of Eq. O |£'i,2) are given 



\D{) = p\D^)+q\D^), 
\D2) - p\D^)-q\D' 



o\ _ „\ nO\ (4) 
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and inversely 



m = ^^{\D{) + \D2))., 
\D') = i;[\Di)-\D^))., 



2q 

where the complex quantities p and q satisfy 

2 



Ml* 



2 2^ 12 



py Mi2-|ri2 



(5) 



(6) 



where M12 and are the complex off-diagonal elements of the 2x2 matrices 
M and F, respectively. In the limit of CP conservation, Di is CP-even and D2 is 
CP-odd0 The eigenvalues of Eq. [3]are: 



71,2 = Ml, 2 — 2^1,2 



Mil - :^rii ± ^ M12 - -ri2 
2 p \ 2 



The eigenstates of Eq. [3] develop in time as follows!!^ 



\Di,2it)) = ei,2W|i?i,2(0)), ei,2W = exp 



-I [ Mi,2 - -ri,2 1 t 



(7) 



(8) 



Using Eq. HI Eq. [S] and Eq. [51 the proper time evolution of a state which is initially 
a pure D'^ is given by: 



where: 



\D\t))^g+{t)\D^) + ^g^{t)\D^), 
P 

\D^{t))=g+{tW) + Pg^{t)\D°). 

q 



g±{t) = [eiWiesW] /2. 



(9) 
(10) 

(11) 



The probabilities for obtaining a D''^ or Z?" at proper time t, starting from an initially 
pure or are: 



|(i?"|i?"(t))|^ = |(i?"|i?"(0)| 

i2 1^-rt 
2 



9+{t)Y 
2 



e ' " |cosh(j/rt) -I- cos(xrt)] , 



(12) 



l5-(^)l' 



[cosh(?/rt) - cos(a;rt)] , (13) 



cosh(?/R) - cos(a;R)] , (14) 



If both X and y are zero, then the probability for a D° to mix to a I?° or for a to 
mix to a will be identically zero for all proper times. If either x or y is non-zero, 
then D^-D'^ mixing will occur. 



=-We use the CP phase convention: CP\D°) = +\D°) and CP\D°) = +\D°). 
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Mi2 and determine the mass and width sphttings AM and AF, respectively: 



AM EE Ml ~ A/2 = 2Re 



AF = Fi - F2 = -41m 



-(M12 - 
P 

-(M12 
P 



(15) 
(16) 



and therefore the characteristics of D^-D'^ mixing. We show the unmixed and mixed 
intensities as a function of the dimensionless variable, Tt, for initially pure states of 
K^, D'^, B° and Bg, in Figs. [Slja-d), respectively. Of the four lowest-lying neutral 
pseudoscalar meson systems, the D'^-D^ system shows the smallest mixing, as noted 
earlier. In the K'^ system, both |x| and \y\ are both of order 1; in the system, 
\x\ and \y\ are both of order 1%; in the and B^ systems, |a;| » \y\. 





rt rt 



Fig. 3. The unmixed (blue) and mixed (red) intensities for an initially pure (a) K*^; (b) Z)"; (c) 
-B"; (d) Bs state. The vertical scale in (b) is logarithmic, the others linear. The values of the mixing 
parameters as defined in Eos. [TlandlSlare obtained using data from Ref. 1191 assuming \ \q/p\ = 1. 



From Eq. [9] (Eq. [TOl), the amplitude that a D° (D") produced at < will 
develop into a linear combination of and and decay into / (/) at time t is: 

{f\n\D\t))^Afg+{t)+Af^^g^{t), (17) 
{m\D\t)) = Ajg+{t) + Aj%^{t), (18) 
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where Af and Af are the and Z?" decay amphtudes to a final state f', Ajr and 
A^ are the and D° decay amplitudes to a final state /: 



Af EE {m\D'), 

Af^{f\n\D"), 

Af^{f\n\D'), 
Af^{f\n\D'). 



(19) 
(20) 
(21) 
(22) 



where "H is the Hamiltonian governing weak decays. Written in terms of the decay 
amplitudes, the general expressions for the time-dependent decay rates r(Z3°(i) — 
/) and T{D°{t) J) are: 



T{D\t) ^ /) = 



-Tt 



|A/ncosh(?/rt) + cos(xn)] 



\Af\'^[cosh{yTt) ~ cos(xR)] 



-2Re ( A}Af-J smh{yTt) + 2Im \A}Af- ] sm{xVt) 



(23) 



nD\t) ^ /) = 



|A.-ncosh(2/rt) + cos(a;ri)] 



\Aj\'^[cosh{yTt) - cos(a;R)] 



-2Re [ A*jAj-j s\TLh{yVt) + 2Im yA*^Aj- j sin(2;rt) 



(24) 



To describe the time dependence for the "wrong-sign" (WS) decay such as — )• 
K+TT- {W -J> if "7r+), we rewrite Eq.[13](Eq.[Ml) in terms of of the Cabibbo-favored 
(CF) amplitude Aj {Aj) and the parameter AJ^ {^f)^ where 



pAj 



(25) 



For decay times t satisfying \xTt\, \yrt\ <^ 1, the decay rates are given by: 
T{D^it) ^ /) = 



Re 



{xj')yTt- 



im(A7Mxrt + ^^i^(ri)2 



T{D'{t) ^ /) = 



|Aj|^ - Re (Aj) yVt - Im (A/) xTt + ^-^-^{Tt)^ 



(26) 



(27) 



Under similar conditions, the time-dependent rates for and decaying to a 
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CP-eigenstate f = f = fcp can be written as 
T{D°{t) ^ fcp) = 

r{D°{t) ^ fcP 



1 - Re (A/,J yVt + Im (Xf,,) xTt + \Xj^, f ^!±l!(n)2 



,-rt 



A 



fcp I 



l-Re(A7^^)yR + Im(A7^^).rt 



JCP 



2 ^2 



(28) 



(29) 



where the terms proportional to e~^*(rt)^ are due to mixing, those proportional to 
e~'"*rt are due to the interference between mixing and decay, while those propor- 
tional to e~^* are due to direct decay. Fig. |4] illustrates the two interfering decay 
paths from an initial to to a final state /. 



A, 



J 



mix 



A, 



Fig. 4. Diagram illustrating the interference between the amplitude for direct decay to / and 
that for W mixing followed by D^^ decay to /. 



1.2. CP violation 

There are three different types of CP-violating effects in meson decays^ 

(1) CP violation [CPV) in mixing; 

(2) CPV in decay, also known as direct CPV] 

(3) CPV in the interference between a direct decay, D° — > /, and a decay involving 
mixing, — > Z?" /. 

CPV in mixing and in the interference between mixing and decay is referred to as 
indirect CPV. 

CPV in mixing occurs when the mixing probability of to is different 
than the mixing probability of D° to D° . As can be seen from Eqs. 1131 and 1141 this 
happens if and only if \q/p\ ^ 1. This type of CP-violating effect depends only on 
the mixing parameters and not the final state of the decay. 

As an example, consider the decay — K'^*^^l~vi. The diagram for this decay 
is illustrated in Fig. [5] Within the SM, the I?° must first mix to followed by 



''For a complete review, see Refs. 1171 1201 1211 
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the direct decay Z?" — > K^*^~^l~Pi. There is no direct decay of to the final state 
K^*^~^l~Di in the SM. Therefore, the time-dependent CP asymmetry: 

, ^ dr/dt{D" -> K(*^+l-Di) - dr/dt{D" ~> K'^*)-l+vi) 

dV/dt{D^ ^ K(*)+l-Di)+dT/dt{D^ ^ Ki*)~l+vi) ^ ' 

is equal to the mode-independent quantity: 

which is used to characterize CPV in mixing. The techniques used to analyze the 
decay £1° ^ K'-*^+1-Di for D"-D° mixing are discussed in Sec. 12.3.41 the results are 
presented in Sec. 13.3.41 



r 




Fig. 5. Diagram for the decay Z?" — > K^*^~^l ui. 



CPV in decay occurs when the amplitude for a decay and its CP conjugate 
process have different magnitudes: ^ 1. In charged D modes where no 

mixing can occur, CPV in decay is characterized by non-zero values for the time- 
integrated asymmetry: 



(33) 



Consider the case where two amplitudes: 

Aj = |ai|e*('^i-'*i) -I- |a2|e^(*^-'^2) 

mediate the decay, where (5i_2 and <f>i^2 are the strong and weak phases, respectively, 
of amplitude yli.2. The weak phase changes sign under CP, whereas the strong phase 
does not. The CP asymmetry Af± can then be written as: 

- 2|ai||a2| sin((?2 - Si)s\n{(j)2 - 0i) 

|aiP + |a2P + 2|ai||a2|cos(52-'^i)cos(02-<^i)' ^ ^ 

Thus, direct CPV will occur {Af± ^ 0) only if the differences between the CP- 
conserving strong phases and the differences between the CP-violating weak phases 
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of the two contributing amplitudes are not zero or multiples of tt. In ne utra l D 
decays, direct CPV is characterized by the mode-dependent parameter AuffS 



Ad ^ " J' ^' . (35) 



\Af/Af 


2 




\Af/Af 


2 
+ 


Af/Af\ 



CPV in the interference between a decay without mixing, D'^ — > /, and a decay with 
mixing, Z?" /, where / can be reached from both and decays, can 

also occur. Consider the time-dependent rate asymmetry of neutral meson decays 
to a final CP eigenstate: 

= r(z^"(t)^/cp)-r(g°(t)^/cp) 

If both CPV in mixing and decay are absent, then jg/pj = 1, \Af^p / Af^p\ = 1 
and therefore, \Xfcp\ = 1- As can be seen by comparing Egs. B51 and Acp[t) 
will be nonzero when AJ^^ ^ ^fcp' Therefore, CPV in the interference between 
mixing and decay will be present if ImA/^p ^ 0, which implies sin (j)f^p ^ 0, where 
Xfcp = |A/cj= I exp(i(/)/pp). The phase (/>/pp is the sum of the phase difference between 
q and p, 

ip = arg (q/p) (37) 

and the (weak) phase difference between Af and Af. In general, the weak phase 
component of arg(A/) is said to characterize CPV in the interference between mix- 
ing and decay. 

The quantities Xj^ and Aj can be evaluated for the modes / = K^n^ and 
/ = K^K^ in terms of Am and the mode-dependent quantities Ru, Ajj, A^^ , 
(j)Kii and (t)KK^ 




(1 


+ Ad)(1 


- Am) 




-Ad){1 


+ Am) 


(1 


-Ad){1 


+ Am) 


(1 


+ Ad)(1 


- Am) 



^ (l-ylD)(l + ^M) ^ ^ 

SK^-<t>K^)^ (39) 



' (1-Ag^)(l + AM 
(1 + Ag^)(l-AM 

For / = K^TT^ the quantity Rd is defined to be 



Rd^^R^Rd, Rt,^\AK^.-/AK+^-\^, R^ = \Ak-^+/Ak-^+\^. (41) 

In the absence of direct CPV, 

= I^A'+.-Mk-.+ P, = |A^-.+ M;^+.-|2, (42) 

are the ratios of the doubly-Cabibbo-suppressed (DCS) to Cabibbo-favored (CF) de- 
cay widths for and decays, respectively. These ratios are O(tan''(0c)) ~ 0.3%, 
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where dc is the Cabibbo angle. AdditionaUy, Sktv is the relative strong phase be- 
tween Aj^+„- and Aj^+^- . The quantities and (j) are in general mode-dependent. 
However, assuming that SM tree-level amplitudes dominate the dec ays, (j!) appearing 
in Eqs. [3i |39] and |40] will be the same for the Ktt and KK modes.12112^ 

Traditional SM estimate s for CP asymmetries in D-meson decays are small, 
less than C'(0.01%).^^^I^3ES1 This is because, to a very good approximation, only 
two generations of quarks are involved in charm mixing and decay, while the 
CKM mechanisrrPfi! requires three quark generations to produce CPV^^ Present 
experim enta l uncertainties on time-integrated CP asymmetries in D decays are 
0(0.1%).!^ Through 2010, all measured CP asy mmetries in D decays were consis- 
tent with zero within experimental errors! I I 

In 2011, the LHCb Collaboration presented evidence for direct CPV by measur- 
ing the difference in time-integrated CP asymmetries between two singly Cabibbo 
suppressed D decay modes: AAcp = Acp{D K^K^) — Acp{D tt+tt^) = 
(-0.82 ± 0.21(stat) ± 0.11(syst))%.'2SI In this difference, the mode-independent in- 
direct contribution cancels. A new standard model calculation of this difference, 
while uncertain to a factor of a few, may accommodate this intriguingly large ex- 
perimental result. New physics, such as supersymme tric gluino-squark loops, could 
also yield direct CP asymmetries as large as 0(1%) 

1.3. Outline 

This paper discusses the experimental status of D^-D^ mixing and CPV as of 
the end of the 2011 calendar year. We review the current results from recent 
coUiding-beam and fixed-target experiments and discuss in some detail the tech- 
niques involved. We survey the primary analysis methods used to study two-body 
and multi-body hadronic and semileptonic decays. Then we present results from 
experimental measurements of mixing and searches for CPV from time- independent 
analyses (those that do not use the proper decay time of the to search for mixing) 
and time-dependent analyses (which do use the proper D'^ decay time) as well as 
quantum-correlated decays. Finally we review future prospects for measurements in 
the near- and longer-term future and summarize the overall status of D^-D^ mixing 
experiments. 

2. Analysis Techniques for Measuring Charm Mixing and CP 
Violation 

2.1. Time-independent Methods 

Time-independent methods provide an important technique for measuring D^-D^ 
mixing and searching for CP violation in charm decays. They also yield information 
on relative strong phases between mixed and direct decays for several different 
hadronic modes of interest to mixing studies. Knowledge of the strong phase 5k-k 
between — > tt^ and — > tt^ allows conversion of the observable y' to 
mixing parameter y (see Sec. 12.2.1]) . 
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As the name implies, these methods do not make use of decay-time information. 
Instead, they count the numbers of and Z?" decays to specific modes when 
the D^-D^ pair has been produced in a quantum-coherent, charge conjugation (C) 
eigenstate. Relative numbers of decay modes of both singly-tagged (ST) events, 
where one 13° or is fully reconstructed, and doubly-tagged (DT) events, where 
both the 1?° and the D° mesons are fully reconstructed, provide information on the 
mixing parameters x and y, the strong phase differences 5i for each decay mode i, 
and their DCS decay rates. This method is especially useful when separating the 
individual D° and decay vertices is difficult, as in the case of non- asymmetric 
energy e^e~ colliders. ^^^^ 

At present, these methods have been performeci^l^^Sl only at the 3.770 GeV 
resonance, but could also be done at higher center-of-mass energies through initial- 
state radiation, if the number of ISR photons can be determined (so that the CP 
state of the coherent D^-D^ pair is known). 

2.1.1. Correlated Decays at 3.770 GeV 

In e+e~ collisions at or above 3.770 GeV that produce a D^-D^ pair, the production 
of the pair may be assumed to proceed through a single virtual photon with J^'-^ ~ 
1^^. At 3.770 GeV, the final state will have C = — 1. At higher energies, additional 
pious and photons may be produced!^ 

e+e- ^ D"^" + ?n(7r") + 7^(7) (43) 

where to, rt > 0. Therefore the D^-D^ pair will be produced with C{D^D'^) = 

— 1"+^. When n = 0, C{D^D^) will be —1. The value of to is not a factor since 
C(7rO) = +1. 

Additionally, if the D^-D'^ pair has relative angular momentum I, then 
C{D°D°) = P{D"D°) = -l' where P is the parity operator. Assuming CP is 
conserved, we can write the wavefunction of the D'^-D'^ state in the center of mass 
system (where the mesons have momentum p and —p, respectively) in terms of ei- 
ther the flavor eigenstates and D'~' or the CP eigenstates Di, D2 (with CP = +1, 

— 1 respectively). If n is even (0, 2, . . . ), the produced D'^D'^ state is 

\D°{p))\D"{-p)) \D"{p))\D°{-P)) - mp))\Di{-p)) \D,mD2{-p)) (44) 
which has C = P = —1. Ifnis odd, the produced state is 

\D"{p))\D°{-p)) + \D°{mD"{-P)) = mp))\Dii-p)) mp))\D2i-p)) m 

with C = P = +l. Therefore at 3.770 GeV when both the D° and D° decay to CP- 
eigenstates, they will have opposite CP. If any same-CP decays occur, the number 
produced will be a measure of the rate of charm mixing. 

To connect the number of like-CP and opposite-CP events to the mixing rate and 
strong phase 6i for a given decay mode i, expressions for time-integrated rates for 
ST or DT events can be calculated from decay amplitudes. Observed rates for one or 
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more modes can be investigated simultaneously, with mixing parameters and strong 
phases obtained from a simultaneous fit to all decay modes under consideration. 

As an example, consider the DT decay to {K^ 7r+, 7r+). From a CP = — 1 
D^D^ coherent state, this rate should be zero in the absence of mixing. A short 
calculation'' yields 



1 



2\ 



where Ak^ = (i^-^+|DO), Ak-. ^ {K-Tr+\D°), tk. = \Ak^/AkA. and 6'^ 
strong phase difference between Ak-k and Aktt'- 



Aktt/Ak 



(46) 
(47) 
s the 

(48) 



where we have incorporated the phase convention used by CLEO in the defi- 



nition of S'j^^. Note that if the mixing rate R 



M = {x^ + 2/^)/2 vanishes, then 
Tic:P=-^\K-Ti+,K-TT+) will vanish. A non-zero rate will be an indication of the 
presence of mixing. This rate can be contrasted with the DT (_fir~7r+, if "•"tt") decay 
rate 



TiCP^-^)(^K-^\K+n-) = \A,,-,, 



\Ak-., 



1 - ri,e-2'^- 



24. cos 2^^, 



1 



1 



1 



(49) 
.(50) 



2 

^ y^)/2. Inclusion of 



Comparison of these rates yields information on Rm = [x 
other DT decay mode pairs permits measurement of x, y, 6i, and branching 
fractions. See Table [TJ 

D° final states used by CLEO-cEH are K^tt^, K+K^, n+n', K^^tt^tt^, K^-k^ , 
K^t:^ , K^r], K^uj, and inclusive semileptonic decays Xe^fe, Xe~T'e. Events contain- 
ing neutral D candidates are selected using two quantities, the beam-constrained 
mass M: 



M : 



(51) 



and the energy difference IS.E = — Eq where Eq is the beam energy, Ejj is the 
sum of energies of the 1?° candidate decay products, and po is the I?" candidate 
momentum. Well-reconstructed candidates will have distributions that peak at 
the D° mass in M and at zero in AE. After mode-dependent cuts on AE are 
imposed, ST yields are obtained by fitting the M distribution and DT yields by 
counting events in a signal region in the two-dimensional M distribution. 

Semileptonic decays are reconstructed inclusively, with only the electron re- 
quired to be identified. Electron identification is performed by use of multivariate 
techniques Decays involving K'^ mesons or neutrinos are reconstructed using a 
missing-mass technique only in DT events. 
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Table 1. Correlated and uncorrelated decay r ates fo r ST and DT events used in 
analysis of coherent D^D^ decays by CLEO-cEDSSl Rates are normalized to the 
branching fraction(s) of reconstructed mode{s) (note that a normalization is used 
where is the branching fraction to mode j when no mixing is present). 
S'_I_(S_) denotes a decay to a CP = +1(— 1) eigenstate; e~ denotes a semileptonic 
decay containing e~ . Rates are given to leading order in x, y, and Rws^ the WS- 
to-RS decay rate ratio. EfEects of CP violation are negligible. Charge-conjugate 
modes are implied. 



ST mode 


Uncorrelated rate 


Correlated rate 


K-TT + 


1 + -Rws 


l + Rws 


s± 


2 


2 


DT mode 


Uncorrelated rate 


Correlated rate 


K-7T+,K-TT+ 


Rws 


Rm 


K-n+,K+T7- 


1 + Rws 


(1 + Rws)'^ - 4r cos (5if^ (r cos <5/f ^ + y) 


K-TT+,S± 


^ + Rws 


1 + Rws ± 2r cos Skt, ± y 


K--K+,e- 


1 


1 — ry cos Sk-h — rx sin Sk-k 


s±,s± 


1 





S+,S- 


2 


4 


S±, e~ 


1 


l±y 



Measurements oi x^, y, , rxsinSK-m and r cos 6 k-k are obtained from the ob- 
served ST and DT yields and external branching fraction measurements using a 
least-squares fit .^^Sl The DT yields provide information on mixing and strong phase 
parameters. Use of ST and DT yields simultaneously provides normalization, so 
that independent measurements of the absolute D^-D^ production rate and the in- 
tegrated luminosity are not required. This method is described fully in Refs. 133 
l39l and HOI including event selection and global fit techniques. Quantum-correlated 
results are presented in Section fS.l.ll 

2.2. Time- dependent Analyses of Two-body Decays 
2.2.1. — > K'^TT^ Wrong-sign Analysis 

In the wrong-sign (WS) decay, — > /C+tt^, the if+7r^ final state may 
be reached either through a direct, doubly Cabibbo-suppressed (DCS) decay, or 
through mixing, — > D'^ , followed by the Cabibbo- favored (CF) right-sign (RS) 
decay, — )■ K^-k~ . Since the two processes involve the same initial- and final 
states and are therefore indistinguishable, interference between the two amplitudes 
will occur. For decays to K^-k~ and decays to A'^7r+, we define the WS 
decay rates relative to the RS decay rates as follows: 

From Eqs. [551 123 ISHl and [31] these are given by: 

Rwsit) = Rd + y'^'iR^m + ""'^ {Tt)\ (53) 
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where the mixing parameters a;'+, y'+ {x' , y' ): 



^ 4^ I ± Am (j>±y' sin 0] , (54) 

1 T Am 



y'^ = \ l-V' T a;' sin 0] . (55) 

V J- T Am 

are the mixing parameters measured in the _D° (-0°) decay modes, cj) is the weak 
phase characterizing CPV in the interference between mixing and decay, and the 
parameters x' and y' are related to the mixing parameters x and y through a rotation 
by the strong phase, 5ktt'- 

x' = X cos 6 K-K + y 5 K-Ki (56) 
y' = y cos 5k-k - x sin 5ktt ■ (57) 

The DCS branching fraction for D'^ and decays is related to the direct CPV 
asymmetry parameter Ajj as follows: 



l±A, 



In the limit of CP conservation {An = Am = = 0), Eg . [53l reduces to: 

Rwsit) = Rd + y'VR^{Tt) + ^l±l^(rt)2. (59) 

The relative WS decay rate allows a determination of x'^, y' and Rd, but not the 
strong phase Sktt- For small mixing parameter values, the main sensitivity to mixing 
comes through the interference term which is linear in y'Tt. Tree diagrams for the 
two amplitudes mediating the D'^ — ^ K^n^ decay are shown in Fig. [51 




Fig. 6. Diagrams illustrating two ways to reach the K~^tt final state from an initial D°. Left: 
direct DCS decay, D° K+tt' . Right: mixing, ^ 15°, followed by CF decay, I)" -s- K+tt'. 



Experiments use the slow pion 7r+ in the strong decay D*^ — > tt+Z? to 
tag the charm flavor of the neutral D at productionjfl The charge of the tt^, 

•^Unless otherwise stated, reference to a given decay mode implies reference to its CP-conjugate 
mode as well. 
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together with the charge of the kaon from the decay of the neutral D allows 
the signal sample to be divided into four categories: two WS decay samples, 
TT+Z?'^, D''* A'+TT^ + C.C., and two much larger right-sign (RS) decay 
control samples, tt+D", D'^ — K~7r^ + c.c A simultaneous fit to the RS 

and WS distributions is performed to determine the direct CF and DCS lifetime and 
the parameters of the decay-time resolution model (from the RS and WS samples) 
and the parameters Rd, a;'^, y' (from the WS sample). The independent variables 
of the fit are rriK-n, the reconstructed Kir invariant mass; Am, the D*~^-D^ mass 
difference; i, the reconstructed decay time, and its measured uncertainty, at- The 
variables mxTr and Am are used to separate signal from background. At BaBAR, 
the vertical height of the beam spot is « 6 /xm. This beam spot information is used 
to constrain the location of the D* vertex, thus substantially improving the deter- 
mination of Am and the reconstructed decay time, t. Fig. [7] shows the projections 
of the rriKTr and Am data an d s ignal and background fit functions from the 2007 
384 fb"^ WS BABAR data set.l^ 




Fig. 7. BABAR, distributions of (a): mKn from WS candidates with 0.1445 < Am < 0.1465 GeV/c^ 
and (b): Am for WS candidates with 1.843 < m^^^ < 1.883 GeV/c^. The projections of the signal- 
and background fits are overlaid, where the random tts background sample peaks in rriKn but 
not in Am; the misreconstructed paired with a tts from a D* peaks in Am but not in m^'ir ; 
the combinatoric background sample peaks neither in m^ir nor in Am. Reprinted figure with 
permission from B. Aubert et al, Phys. Rev. Lett. 98, 211802 (2007). Copyright 2007 by the 
American Physical Society.Q2] 



2.2.2. Lifetime Ratio Analysis 

In this section we describe analysis techniques which measure the decay-time distri- 
butions of neutral D mesons decaying to CP eigenstates and CP mixed states. The 
potential of this method was first described in Ref. HU and the first experimental 
results were presented utilizing these techniques in Ref. 1421 In the last ten years. 
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several experimental collaborations have measured D^-D^ mixing and CP violation 
observables with increasing precision by comparing the rate for mesons decaying 
to flavor-specific final states. 

In particular, the mixing parameter y (Eq. [5]) may be measured by comparing 
the rate of decays to CP eigenstates with decays to non-CP eigenstates. If decays 
to CP eigenstates have a shorter effective lifetime than those decaying to non-CP 
eigenstates, then y is positive. 

Experimentally, one is interested in collecting high purity samples with large 
statistics of decays to final states of specific CP content. The two-body SM 
processes with these characteristics are the singly-Cabibbo-suppressed decays to 
CP-even eigenstates, D'^ — > K^K~ and — > tt+tt^, shown in Figs.[5I^a) and[5{b) 
respectively, and the Cabibbo-favored decay to CP- mixed final state, — > i^^7r+, 
shown in Fig. [21 and the corresponding CP-conjugate decay processes. 

Similarly to the two-body final states, the mixing parameter y can also be mea- 
sured by analyzing the CP-odd component of — K^K'^K~ decays, by means of 
comparing the mean decay times for different regions of the three-body phase space 
distribution of the final stateH 




u * u u * u 



Fig. 8. Feynman diagrams for singly-Cabibbo-suppressed decays of D" to CP eigenstates, 
D° K+ K- (a), and -s- vr+vr- (b). 




Fig. 9. Feynman diagram of the Cabibbo-favored decay D —> K n 



"^Details of this measurement will be discussed in Sec. 12.31 
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Neglecting the quadratic (mixing) terms in Eqs. [55] and 1291 an approximation 
valid when |A/pp| ss 1 and \xrt\, \ yTt\ ^ 1, we obtain the the following expressions 
for the time-dependent decay rates T{D'^{t) h+h^) and T{W{t) h+h^): 



) = 
) = 



lA 



h + h 



[Re{Xh+h- )y - ^MK+h- )x\Tt} , 



[Re(A, 



(60) 



T{D°{t)^ K+TT-) = 



-rt 



K+TT- 



(61) 



where h — tt,K. In the absence of direct CP violation (as expected in the SM), but 
allowing for a small indirect CP violation (with a weak phase \(f)\ <Si 1), we can write 
Xf^+f^- = \q/p\ e"^ . To a good approximation, these decay-time distributions can be 
treated as exponentials with effective lifetimes given by Ref. |43] 



'hh 



'hh 



h+h-) 
h+h-) 



= t{W - 

= TK-Tr+ 
= TK-TT+ 



1 + 


q 




p 


1 + 


P 




9 



{y cos (j) 
{y cos 4> 



X sm ( 



X sm ( 



(62) 



as before h = n, K . Combining these quantities we then define the parameters ycPy 
Ar and AY as: 

{Thh) ' 

-jP" h-h+) - t{D" h+h-) 



ycp 

Ar 



-(£>" ^ h-h+) + t{D" h+h-) ■ 



AY 



(Thh) 



At 



- [j^y^ '^^hhll'^ ^^'^ h = Tr,K. In the 
0. In the absence of D^-D^ mixing, 



where (. . .) implies average over fiavors, {rhh) 
limit of CP conservation, ycp = y and Ay : 
both ycp and Ay are zero. 

Measurements of ycp have been conducted at e+e~ colliders [BaBAR, Belle 
and CLEO) as well as fixed-target experiments (FOCUS and E791). Historically, 
experiments at e+e~ colliders have relied on the kinematic separation of charm 
decays at high center-of-mass momentum (from e+e- — > cc) to reduce backgrounds. 
In addition, excellent particle identification and tracking capabilities for hadrons 
over a large range of momenta are required when measuring ycp- 

The BaBAR and Belle experiments have both produced measurements of ycp 
and Ay^"^ 44 45|^y means of selecting highly pure samples with high statistics of 
candidates decaying to K+R- , tt+tt" and K~7r+ final states, as shown in Fig. [101 
In these experiments, D mesons are produced from cc initial states and as sec- 
ondaries from B decays, those produced from cc events are used in the lifetime 
ratio measurements by choosing high momentum D mesons as well as minimizing 
other backgrounds. However, slightly different strategies were followed by different 
experiments, but the description of what follows is in general correct for all ycp 
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measurements. We focus primarily on the measurements done at the _B-factories as 
those are the most precise. 

The so-called tagged technique is implemented by reconstructing D*^ D^T^t 
decays, where the slow tt^ determines the flavor of the decaying neutral meson as 
or at creation. The D° candidates are selected by kinematically combining pairs 
of oppositely-charged and tt^ tracks that have a common vertex, and have an 
invariant mass typically in the range between 1.8 GeV and 1.92 GeV (approximately 
±60 MeV around the nominal Z?" masJ^^. The kinematic flt provides the 1?° decay 
position and its momentum vector pjjo, which is required to point back to the e+e" 
interaction region. The candidate and the slow tt^ are also required to form a 
common vertex in the interaction region. For each candidate the proper decay 
time t and its error at are calculated using the decay length / = /S^ct = ctpuo /ttido. 

To further suppress background events, Belle exploits the distribution of the 
energy released in the D*^ decay given by Q = ™_d* ~ m^o — m.^\ equivalently, 
BaBAR uses the distribution of the mass difference Am of the reconstructed D*^ 
and the candidates in the event. Typically candidates are required to be 
within ±0.1 MeV of the peak of the Am or Q distributions. 

The BaBAR invariant mass distribution of tagged events of different decay 
channels are shown in Fig.lloi^The shaded area shows the sample of events used in 
the lifetime measurements. These events were selected after particle identification, 
tracking, and vertex probability requirements. 

In the tagged sample, the charge of the reconstructed D*^ allows the determi- 
nation of the lifetime separately for or D'^ decays. The lifetimes are determined 
by performing a simultaneous maximum likelihood flt to the reconstructed decay 
time and its error to all five decay samples {D'^ — > A' ^7r+ and — > K^n~ decay 
samples are combined into one). In general, there are three main PDF components 
entering the lifetime fit: signal, combinatoric background, and misreconstructed 
charm decays. 

In the so-called untagged technique, there is no reconstructed D*^ , hence it is 
not possible to identify the initial flavor of the decaying meson. In the untagged 
analysis, the D'^ may be produced directly from a cc state or as a decay product 
of a higher mass resonance (other than D*^). The momentum is required to 
point back to the beam spot in order to reduce backgrounds. In order to exclude 

mesons coming from B decays, D° candidates with momentum in the e'^e~ 
center-of-mass (CM) frame less than 2.5 GeV/c. 

In the untagged BaBAR analysis,!^ all events appearing in the tagged data sam- 
ple are removed from the untagged sample in order to treat the tagged and untagged 
results as statistically independent from one another. The signal yields in the un- 
tagged data samples are about 3.5 times larger than those in the respective tagged 
samples; however, their purity is lower and the systematic uncertainties due to the 
higher backgrounds are more challenging. 

The signal PDF is generally described by an exponential convolved with a resolu- 
tion function, which is composed of three Gaussian functions sharing some common 
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parameters between them. The high statistics of the Cabibbo-favored D'^ — > K^tt^ 
sample drives the determination of the resolution function parameters in the life- 
time fit. In the tagged analysis, the random combinatoric background in the signal 
region is determined from a sideband region in D'^ invariant mass (m^jo) and Am. 
In the untagged analysis, two sideband regions in m^jo are defined, one above the 

mass peak and one below. A small background component corresponding to mis- 
reconstructed charm decays that have long lifetimes and can thus mimic the decay 
time of signal events is included. The proper time distribution for this background 
is taken from Monte Carlo (MC). 

While systematic uncertainties are expected to cancel in the lifetime ratio, the 
sources of backgrounds are different for each final state, hence systematics from 
background sources are not necessarily expected to cancel. The main signal model 
systematic uncertainties are the selection of the invariant mass signal region 
window (central position and size), opening angle distributions, and variations of the 
signal resolution model. The systematic uncertainties associated with backgrounds 
are the combinatorial PDF model and its normalization, and the misreconstructed 
charm PDF model (taken from simulation) and its normalization. 

Results from these measurements are discussed in section [3.2.21 below. 

2.3. Time- dependent Analyses of Hadronic Multi-body Decay 
Modes 

Amplitude analyses of multi-body D'^ decay modes provide what are potentially 
the most definitive measurements of charm mixing parameters. Advantages include 
the ability, for some decay modes, to measure mixing without the ambiguity of an 
unknown strong phase or insensitivity to the sign of x' that limits the measurement 
to x' and y' rather than x and y, as is the case with the time-dependent analysis 
of K^TT^ decays. Multi-body decays useful in this regard include ^■ 

K^K~^K^ or iC^TT+TT^, which we will generically designate as K'^h^h" where h 
represents K or tt. Three-body decays also include Z?" — X+tt^tt". Four-body 
decays include K^tt^tt^tt^ . Three-body decays are amenable to "Dalitz- 

plot analysis," while higher-order decays require other methods. 

2.3.1. D° K°h+h- Analysis 

BaBAR, Belle, and CLEO have performed studies of D^-D^ mixing using the decays 
ZJO ^ A'Ott+tt-, do K^^K+K-, or both.46 47 48 r^j^g -^^^^ jg ^-j^g Dalitz- 

plot distribution of selected decays using the time-dependent formalism given 
in Eq. [23] (for £>°) and Eq. [24] (for D°). The variation of the decay amplitudes 
Af, A J and their conjugates across the Dalitz plot must be taken into account. 
We define A(s+,s„) to be the amplitude for — >• Kgh~^h~ and A{s+,S-) to 
be the amplitude for D'^ — > Kgh^h^, where s+ and s_ are the coordinates of 
a given position in the Dalitz plot, e.g., s+, s_ = to^o^_ , 7n^_^_|_ (CLEO) or s+. 
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Fig. 10. The BABMl, reconstructed D'^ mass distributions for the three tagged samples, within 
±0.8 MeV/c^ of the peak of Am. The shaded region indicates the events used in the Ufetime fit. 
(The structures appearing above 1.92 GeV/c^ in the K' K+ decay mode, and below 1.81 GeV/c^ 
in the 7r+7r~ decay mode, are mainly due to candidates with misidentified kaons or pions.) Also 
shown are the yield and purity of the three samples as calculated inside the ±15 MeV/c^ mass 
window used in the lifetime measurements. Reprinted figure with permission from B. A uber t et 
al, Phys. Rev. D 78, 011105{R) (2008). Copyright 2008 by the American Physical Society.El 



S- = TO^o^+, m^o^- (Belle, BaBAR). In order to fit the Dalitz-plot distribution as a 
function of time, it is necessary to assume a Dalitz fit model. These models typically 
include a coherent sum of ten to twelve quasi-two-body intermediate resonances 
plus a non- resonant component. P- and Z?-wave amplitudes are modeled by Breit- 
Wigner or Gounaris-Sakurai functional forms, including Blatt-Weisskopf centrifugal 
barrier factors. In the i3ABAR analysis, the S'-wave dynamics are modeled using a K- 
matrix formalism (tttt), a Breit-Wigner plus non-resonant contribution (-R'tt), and 
a coupled-channel Breit-Wigner model describing the ao(980) and Breit-Wigner 
models for the /o(1370) and ao(1450) [KK). 

As far as mixing is concerned, the interesting Dalitz plot regions are where the 
CF and DCS amplitudes interfere and regions where CP eigenstates predominate. 

These analyses proceed in a manner similar to the two-body, time-dependent 
analyses: they make use of the sign of the slow pion tt^ from a D*"*" decay to 
tag the neutral meson as or 13° at its creation. After selecting appropriate- 
quality charged tracks, tt+tt^ pairs that have an invariant mass close to the 
mass (typically within 10 MeV) are selected, forming a candidate. Another set 
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Fig. 11. BABAR. distributions of m^jo and Am for the Kgh^h~ analysis, after all selection cuts. 
The shaded scale indicates the number density of events. The boxes indicate the signal regions as 
used in the analysis for the mixing fits, (a) K^tttt. (b) K^KK. Reprinted figures with permission 
from P. del Amo Sanchez et al, Phys. Rev. Lett. 105, 081803 (2010). Copyright 2010 by the 
American Physical Society^S' 



of oppositely-charged tracks that share a common vertex are combined with the 
candidate to form a Z?" candidate. This ahows for the K'^ decay vertex to be 
displaced from the decay vertex. A kinematic fit then provides the decay 
position and its momentum vector Pdo, which is required to point back the the 
luminous interaction region. The decay time t is calculated using the decay 
length I — P^ct = ctp£)0 /niDa, along with its error at, on an event-by-event basis. 

Background sources include the random tt^ background, where an incorrect as- 
signment between a low-momentum pion and a good decay has been made, 
misreconstructed -D", and combinatoric background. A few other sources of back- 
grounds (classification varies from experiment to experiment) may also be included 
in the fit model as well to model specific non-signal decay modes. 

The time-dependent analysis uses candidates from a two-dimensional signal re- 
gion of Mj^o, the reconstructed candidate mass, and either Am (BaBAR), or 
Q ~ rrij^o^^^^ — TOj^Oj,.^ — mTT, the available kinetic energy released in the £)*+ 
decay (Belle). See Fig. [TT] BaBAR (Belle) determines the yields of signal and back- 
ground in the signal box by fitting the mjyo and Am (Q) to PDFs characteriz- 
ing each background source over the full range in m/30 and Am {BaBAR) or Q 
(Belle), and rescaling the component yields to the signal region. Belle finds 534,410 
signal candidates in 540 fb"\SSI BaBAR finds a signal yield of 540,800 (79,900) 
K^TTTT (K^K+K-) with purity 98.5% (99.2%) K^^rnr (K^^K+R-) in 468.5 fb~^ 
of data.l^ Fig. [T^ shows the Belle experiment's time-integrated distribution of 
decays and projections of the fit to the data where mj_ ~ rn'^j^o^^ for D'^ decays 



September 27, 2012 0:11 WSPC/INSTRUCTION FILE D' mixing 'review 



22 Chavez, Cowan, and Lockman 



It 

o 

CM 
> 

O 

CM + 




m? (GeV^/cl 



ml (GeV^/cl 




(GeV^/cl 



m 2 (GeV^/c^)' 



Fig. 12. Dalitz plot and fit projections for the Belle — )■ K^n~^n analysis. In the projections, 

for decays, and 



the data are shown as points and the fit as a solid line. mj_ = m'^Q^j_ 
7Tij_ = m'^g for D" decays. See Ref. 1481 for details of the 18 quasi-two-body resonance plus 

non-resonant background Dalitz model and resulting fit. Reprinted figure with permission from L. 
M. Zhang et ai, Phys. Rev. Lett. 99, 131803 (2007). Copyright 2007 by the American Physical 
Society. 



and mj_ — m^o^^ for D'^ decays SHI, 

To determine the mixing parameters, PDFs are defined that inchide the depen- 
dence of the signal and background components on decay time t, at, and location 
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in the Dalitz plot. Included in the signal PDF are the matrix elements 

M{s.,s+,t) ^ A{s^,s+)g+{t) + ^Ais-,s+)g-it), (63) 

M{s^,s+,t) = A{s_,s+)g+{t) + ^A{s^,s+)g_{t), (64) 

q 

which are convolved with a decay-time resolution function that depends on position 
in the Dalitz plot. Eqs . l63l and l64l are generalizations of Eqs. [IT] and [18] to multi-body 
decays. Different resolution functions are used for K^tttt and K^KK distributions. 
The decay-time resolution function is a sum of Gaussians with widths that may scale 
with the event-by-event, decay-time error 5t, and also depends weakly on position in 
the Dalitz plot. Belle uses three Gaussians with different scale factors and a common 
mean, which are allowed to vary in the fit. BaBAR uses two Gaussians that scale 
with St, one of which is allowed to have a non-zero mean (to offset), and a third 
Gaussian which does not scale with St. The results of this procedure are discussed 
in section [3. 3. II 

2.3.2. D" K+tt-tt" Analysis 

As in the case of the two-body WS decay Z?" — K^tt~ , the three-body WS decay 
— > iiT+vr^Tr" can occur through DCS decay or via — >■ D° mixing followed by 
the CF decay D" K+n-jr". With a WS branching fraction of (3.04±0.17) x 10"'* 
compared with (1.47 ± 0.07) x lO""* for D° K-tt+^ the Ktttt^ channel is 
competitive in sensitivity to the two-body channel, despite the lower efficiency of 
reconstructing the three-body final state. 

Reconstruction details of — Ktttt'^ events vary from experiment to experi- 
ment, but the basic selection process is as follows. The decay — > tTsD" is used 
to tag the flavor of the neutral D at production. To form Z)" candidates, pairs of 
oppositely charged tracks originating from a common vertex are combined with a 
7r° candidate whose momentum in the laboratory is > 300 MeV/c. D° candidates 
whose mass is within ~ 60 MeV/c^ of the nominal mass are retained. Particle 
identification requirements are imposed to reduce feedthrough of doubly misiden- 
tified, CF candidates into the WS sample. The momentum of each candidate 
is required to point back to the interaction region and its momentum in the CM 
system (p*) is required to satisfy p* > 2.5 GeV/c to suppress candidates from 
B decay. 

Each I?" candidate is paired with a slow pion tTs to form a D*~^ candidate. 
D*~^ candidates which have an appropriate value of Am (or, equivalently, Q; see 
Section [2. 3. 1|) and have sufficiently good per degree of freedom from the kinematic 
and/or vertex fits are retained. 

Background sources considered are random tTs (an incorrectly associated tTs com- 
bined with good forming a D*^ candidate), incorrectly reconstructed charm 
decays, and uds combinatorial background. Maximum likelihood fits to the two- 
dimensional {mxTrTrOy Am) distribution are performed to determine the yields of 
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Fig. 13. Belle Ktttt^ reconstructed mass (nij^^^a) and available kinetic energy (Q) distributions, 
showing signal and background components as determined by a fit to the two-dimensional (mj^^^o , 
Q) distribution. Top row: RS; bottom row: WS. Reprinted figure with permission from X.C. Tia n 
et al, Phys. Rev. Lett. 95, 231801 (2005). Copyright 2005 by the American Physical Society.HSl 



signal and background candidates in the RS and WS samples. As an example, see 
Fig. m for the Belle fit to the D° candidate mass and Q distributions. 

BaBAR analyzed the D° — >■ Ktttt'^ mode using two different methods. The first, 
method I,^uses the different decay-time dependence of DCS and mixed decays and 
analyzes regions of phase space chosen to optimize sensitivity to mixing. Although 
the rates of DCS and CF decays vary across the Dalitz plot, the mixing rate is the 
same at all phase space points. 

The time dependence of the WS-to-RS decay rate ratio can be expressed for a 
given phase-space region (a tilde indicates integration of a quantity over this region) 
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^^1^ = + ay'./R^iTt) + ^^^(^^)^ (65) 

^ RS \'') ^ 

where a is an averaging factor that accounts for the variation of the strong phase 
over the phase space region (0 < a < 1). Rd is the DCS branching ratio, x' and y' 
are the mixing parameters x and y rotated by an integrated strong phase S: 

x' = a;cosi5 + v sin(5, 

-/ • ? ? (do) 

y = —X sm d + y cos o . 

Note that Rm — (i'^ + y'^)/2 = (x^ + y^)/2 is independent of the integration region. 

After signal and background yields are determined, a fit is performed to the 
decay-time distribution. The functional forms of the PDFs are determined based on 
MC studies, but all parameters are determined by fitting the data. The large RS 
signal is used to determine the resolution function used in the WS fit. The observed 
decay-time dependence of the WS PDF is given by Eq. 1^5] convolved with the 
resolution function. Projections of the reconstructed Z?*^ mass, Am, and decay-time 
fits are shown in Fig. [TH 

The second method (method 11) used by BARAfPSl to search for mixing in 1)° — >■ 
K-K-K^ decays is a time-dependent, Dalitz-plot analysis that uses an isobar model^ 
to describe the dynamics. The time-dependent decay rate for WS decays to a partic- 
ular final state / at a given point in the Dalitz plot (si2,si3) = ("^/f+^- j '™/f+^o)j 
and assuming |y| ^ 1, may be written as 

d7V/(Si2,Si3,t) _ 



dsi2dsi^dt 



e 



13t 

\Af\'' + \Af\\Af\{ycos6f -xsin5f)Tt + ^l±^|A^|2(rt)2 



(67) 



where / = K+tt-tt^, the DCS amplitude is A/(si2,si3) = (/|H|£i°), the CF am- 
plitude is A/(si2,si3) = {f\n\W), and (5/(si2,si3) = m:g{A*/Af). 

Written in terms of normalized mixing parameters and normalized amplitude 
distributions, the time dependence can be expressed as: 

dNf{si2, si3,i) 



oc 



dsi2dsi3dt 

where A^'^^ = ^//y / \Af\'^dsi2dsiz, = Aj/J J Ajl'^ds^dsis are normalized 



distributions, rp = y / \Af\'^dsi2dsi3/ J \ Af\'^dsi2dsi:i, and x = x/vq and y — y/rQ 
are normalized mixing parameters. 

The isobar model parametrizes the amplitudes A f and as a coherent su m of 
seven resonances plus a Kn S-wave component derived from Kn scattering data,^^ 
including a non-resonant component. The high-statistics RS sample (~ 659, 000 
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Fig. 14. BABAR fit projections for Kmr^ analysis method I (see text), showing data (points 
with error bars) and PDFs (dotted fines), (a) Reconstructed mj^^^o for candidates satisfying 
0.1444 < Am < 0.1464 GeV/c^; (b) Am for candidates satisfying 1.85 < m^^^o < 1.88 GeV/c^; 
(c) decay-time t satisfying both mass selections in fa) an d (b); and (d) signal-enhanced version 
of (c) using a channel-likelihood signal projectionl5lIS2l Reprinted figure with permission from 
B. Aubert et al, Phys. Rev. Lett. 97, 221803 (2006). Copyright 2006 by the American Physical 
SocietylSE 



candidates) is used to determine the isobar model parameters for CF decays and 
the decay-time resolution function for both the RS sample and the WS sample (~ 
3000 candidates) . See Fig. 1151 Sensitivity to the mixing parameters arises primarily 
from the interference terms (linear in t) in Eq. [67] and Eq. 1681 PDFs expressing 
the dependence of the WS decay rate on Dalitz plot position and decay time are 
convolved with the decay-time resolution and a fit performed, determining the DCS 
isobar model parameters (amplitudes and phases) and the mixing parameters. 
An unknown strong phase difference ^KTrir" between the DCS decay — > 
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mj^.^.^o [GeV/c^] A m [GeV/c^] 

Fig. 15. BABAR RS (top left) and WS (top right) Dalitz distributions and reconstructed WS 
mass (bottom left) and Am (bottom right) distributions for — >■ Kmr^ using method II. 
Shaded regions indicate signal (white), mistag background (gray), and combinatoric background 
(dark gray). Reprinted figure with permission from B. Au ber t et al, Phys. Rev. Lett. 103, 211801 
(2009). Copyright 2009 by the American Physical Society.EHl 



and the CF decay D'^ — > p cannot be determined in this analysis, so the mixing 
parameters measured are 

Results of the two methods are discussed in Section 13.3.21 
2.3.3. D" K+TT-TT+TT- Analysis 

The D'^ — > K^TT^TT^TT^ (CF K3tt) decay has been used to study charm physics 
since soon after the discovery of the Z?+ and mesons. An early search fo r w rong- 
sign decays saw no significant signal, but did set limits on the WS rate.l^E791 
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reportecl^ a measurement attributed to DCS decay of R^^"" = [0.25t°;34 ± 0.03]% 
by analyzing the distribution of wrong -sig n D° decay times. In a time-integrated 
measurement, CLEO reported evidence^ for wrong-sign _D" if+7r^7r+7r~ de- 
cays. They found a 3.9 standard deviation result in the wrong-sign to right-sign 
branching fraction: R^^^ = [0.4lt[! :[^(stat) ± 0.04(syst) ± 0.10(phase space)]%. 

The Z?" — > K3tt decay offers some advantages over Ktt, and a couple of 

difficulties. One advantage is that the RS branching fraction for Kmnr of 

~ 8.1% is twice that of D'^ -> Ktt of ~ 3.9%. Another is that the vertex resolution 
of the four-body K3n decay is usually better than that of the two-body Kn decay, 
which leads to an improved decay-time resolution. These advantages are somewhat 
offset by the reduced efficiency of reconstructing the four-body decay relative to the 
two-body decay and by complications in determining the mixing parameters x and 
y due to variations in the strong phase over the four-body phase space (the mixing 
rate however, is independent of position in phase space). As in the case of 

— Kmr'^, the strong phase SKTnrTr cannot be determined in this analysis alone. 
To date, no amplitude analysis of — > K3tt decays has been attempted. 

2.3.4. Analysis of Wrong-sign Semileptonic Decays 

The WS semileptonic decays D'^ ~> K^^^^e^v^ and D'^ — > K^*^^^^^^ offer unique 
features to searches for D^-D'^ mixing. One unique feature is that doubly Cabibbo- 
suppressed, wrong-sign decays do not occur in the semileptonic mode in the SM. 
This simplifies the time-dependence of the WS rate relative to the RS rate, Eq. EH 
to: 

Rws{t) = ^{Ttf. (70) 

The WS decay rate is thus directly sensitive to the presence of mixing, as there is no 
contribution from either DCS decay or from interference between DCS decay and 
mixing. On the other hand, semileptonic decays present a challenge not encountered 
when analyzing hadronic decays: the presence of the unobserved neutrino in the final 
state precludes exact determination of the D° candidate mass, leading to degraded 
decay-time and mass-difference resolutions and higher backgrounds. 

Distinguishing characteristics of mixing in WS semileptonic decays include the 
quadratic time dependence of Eq. [70] and a peak in the available kinetic energy 
spectrum Q = m(KlvKs) — m{D'^) — raii:) near 5.8 MeV/c^ (or, equivalently, in 
Am = m{KlvKs) — m{Klv) near 145 MeV/c^). WS semileptonic decays share the 
peaking behavior in Q (or Am) with RS semileptonic decays, but have a time 
dependence modified by the quadratic term given in Eq. [70] instead of the pure ex- 
ponential decay-time distribution characterizing the RS decay. Semileptonic decays 
are also susceptible to feed-through from RS — )■ Ktt decays, where the kaon is 
mis-identified as a lepton and the pion as a kaon. This is particularly a concern 
in the case of semi-muonic decays, since the kaons and pions are more prone to 
mis-identification as muons than as electrons. 
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Many searches for mixing using semileptonic decays have been carried out. The 
K^u mode in particular has been used to sear ch for and set limits on mixing since 
shortly after the discovery of the meson. I -^^I^ More re cently, E791, CLEO, 
Belle, and BaBAR have reported measurements. l^lESIMEllM! 

The E791 analysis estimates the missing momentum of the neutrino by using 
the measured decay vertex positions of the D*~^ and the the kaon and lepton 
momenta, and attributing the D° mass to the secondary decay. This results in a two- 
fold ambiguity, which is resolved by always choosing the higher-momentum solution 
for the (motivated by MC studies). This choice results in some degradation in 
the decay-time resolution which is accounted for as a systematic error. 

The Belle analysis applies the following procedure to estimate the neutrino four- 
momentum and consequently. Am. Applying four momentum balance to the initial 
e+e~ system, the Kl system, the missing v, and the rest of the event, an approxima- 
tion for the missing momentum is obtained. This value is refined by the use of two 
additional constraints. First, a D*^ mass constraint is applied to the Klv system, 
resulting in a scale factor that is used to produce a refined m{Klv'Ks) value, the 
ra{Klv) mass having been fixed to ra{D'^). A second constraint on ra^ is applied, 
resulting in a correction to the angle between the three-momentum of the Kl system 
and that of the rest of the event. 

BaBAR has published two semileptonic mixing analyses, one using a single-tag 
(£)*+ D'^irf) method and the other using a double-tag method. The single-tag 
analysis includes both Kev and — > K*ev decays, and treats them essen- 

tially the same way. No attempt is made to reconstruct the K* explicitly; its kaon 
daughter is used directly in reconstructing the 1?°, as if it were a 13° daughter. After 
selection cuts are imposed, resulting K, e, and tt^ tracks, the position of the K-e 
vertex, and the event thrust axis are used to reconstruct the three components of 
the momentum vector by means of three neural net estimators. These estima- 
tors have been trained using 0(10^) simulated signal events to reproduce the 
momentum vector components. Events are required to pass a neural net selection 
which discriminates prompt charm from background events. The majority of the 
remaining background comes from charm events not from B decays where a ran- 
dom charged pion has been combined with a charged K daughter and an electron 
daughter from the charm decay, or with K and electron combinations not from 
a common parent. Understanding the origin of these backgrounds is important as 
they do not share exactly the same decay-time distribution as true charm decays, 
and this must be accomodated in the decay-time fit. After performing an extended 
maximum likelihood fit to the large RS data sample, which determines many of the 
PDF parameters describing the RS and WS Am and decay-time PDFs, the mix- 
ing quantities are determined from a fit to the WS data, including the decay-time 
information. 

The BaBAR double-tagged analysis attempts to address the predominant back- 
ground present in the singly-tagged analysis: feedthrough into the WS sample from 
RS semileptonic decays, where the has been wrongly associated with a ran- 



September 27, 2012 0:11 WSPC/INSTRUCTION FILE D' mixing 'review 



30 Chavez, Cowan, and Lockman 



dom slow pion tTs. In addition to the D*'^ tag, a second tag is constructed by 
requiring either a fully reconstructed, high- momentum, hadronically decaying D'^ 
or in the the opposite hemisphere. While greatly improving the purity of the 
tagged sample, the selection efficiency drops by more than a factor of ten. Addi- 
tional background suppression criteria are imposed which bring the sensitivity of 
this analysis to about the same as the single-tag analysis above. Resulting Am 
distributions for RS and WS events are shown in Fig. [T6l 




Fig. 16. BABAR semileptonic Am distributions from Ref . 1641 showing both singly-tagged (white) 
and doubly-tagged distributions (gray) Left plot: RS data (points) before the double-tag selection, 
along with the total fit projection (solid line), and the background fit projection (dashed line). Inset: 
RS Am distribution after applying the double-tag selection. Right plot: WS data events satisfying 
all selection criteria (gray histogram) and all but the double-tag selection (white histogram). 
Reprinted figures with permission from B. Aubert et al., Phys. Rev. D 76, 014018 (2007). Copyright 
2007 by the American Physical Society. 



3. Current Experimental Results 
3.1. Time-independent Experiments 
3.1.1. Correlated decay results at 3.770 GeV 

Using methods described in Section 12.1.11 the CLEO Coll aborati on reported the 
first measurement of the strong phase difference Sktt in 2008.1^1^^. From 281 pb^^ 
of data collected at ^/s = 3.770 GeV with the CLEO-c detector, the correlated anal- 
ysis was performed using different sets of external measurements as input. These 
included: measurements of two-body branching fractions; the previous, plus 
measurements of the time-integrated WS rate Rws = r(i)" K^-k^) /T{D^ ^■ 
K^t:^) and the mixing rate Rm (the "standard" fit); and the previous, plus mea- 
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surements of x, y, x'^ , y' , and (the "extended" fit). Correlations between external 
inputs were incorporated in the fits. In the standard fit, xsiiiSktt is fixed to zero. 
In the extended fit, this condition is relaxed. 



3950108-001 




-0.5 0.5 1.5 0.0 0.5 1.0 1.5 2.0 



cos5 

Fig. 17. Fit likelihood from the CLEO standard fit, including statistical and systematic uncer- 
tainties, (a) For cos 5; the hatched region contains 95% of the physically allowed area, (b) Two- 
dimensional likelihood for cos<5 and y. Reprinted figure with permission from J. Ros ner et al., 
Phys. Rev. Lett. 100, 211801 (2008). Copyright 2008 by the American Physical Society.EH 



Systematic uncertainties accounted for in the analyses included estimates of ef- 
ficiencies for particle identification, reconstructing tracks, and for reconstructing 
neutral neutral decays (if^ and tt"). Other systematic uncertainties included ef- 
ficiencies for T] reconstruction, selection cuts, fit model description, and detector 
and physics modeling. The largest systematics were for rj reconstruction (4.0%) and 
AE = Ejjo — Ehcam selection (0.5-5.0%). Bias estimates on the fitting procedure 
were obtained by studying a sample of simulated D^-D^ decays fifteen times the size 
of the recorded dataset. Biases from the fitting procedure were less than one- half 
the size of the statistical errors on the fitted parameters. 

Results from the standard and extended fits are given in Table [2l Likelihoods 
from the standard fit are shown in Fig. [17] and from the extended fit are shown in 
Fig. 1181 The final result for Sx-m including asymmetric errors estimated from the 
shape of the likelihood function shown in Fig. [iTl is cos 5k-k = l.O^tol) ± 0.06. 
Limiting to the region |coS(5/<-7r| < 1, \5kti\ < 75° at the 95% confidence level. 
From the extended fit, they obtain cosSktt = 1-10 ± 0.35 ± 0.07 and xsinSK-rr = 
(4.4^" ± 2.9) X 10-3 and Sktt = (22i^2iii)°- In both cases, the statistical errors 
were obtained by inspection of the log likelihood. 

In a recent preliminary analysis, - CLEO extended its quantum correlated co- 
herent decay analysis to measure the strong phase differences in D*^ — > JC+vr", 
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Table 2. CLEO mixing and 5i^^ strong phase difference measurements. 
From Refs. |3I1 and I35l See text for fit descriptions. 



Parameter Standard fit Extended fit 



y 


(-45 ±59 ±15) X 10-3 


(6.5 ±0.2 ±2.1) X 10-3 




(-1.5 ±3.6 ±4.2) X 10-3 


(0.06 ±0.01 ±0.05) X 10-3 


r2 


(8.0 ± 6.8 ± 1.9) X 10-3 


(3.44 ± 0.01 ± 0.09) X 10-3 


cos Sk-k 


1.03j;o:3J. ± 0.06 


1.10 ±0.35 ±0.07 


X sin &K-K 


Fixed at 


(i.4+];l ±2.9) X 10-3 






/Qr, + ll+9 \0 
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Fig. 18. Fit likelihoods from the CLEO extended fits, including statistical and systematic uncer- 
tainties. Hatched regions contain 95% of the physically allowed area, (a) costS. (b) xsiniS. (c) S. 
(d) Two dimensional likelihood for cos<5 and xsin5. Reprinted figure with permission from J. 
Rosner rf al, Phys. Rev. Lett. 100, 211801 (2008). Copyright 2008 by the American Physical 
Society. -^Ij 



D" if+TT-Tr", £»" ^ K+TT-TT+TT- and D° K^j^h+h- , h = Tr,K, using the 
full dataset (818pb^^) together with additional single- and double-tag modes. This 
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analysis makes direct measurements of and sinSKTr, resulting in approximately 
a factor of two smaller (and more symmetric) statistical uncertainties on cos Sk-k- In 
the near future, BES-III will likely produce strong phase difference measurements 
with improved statistical precision. 



3.2. Results from Time- dependent Analyses of Two-body Decays 
3.2.1. D'^ — > K^TT^ Wrong-sign Decay Results 

Several experiments, E691,ISSI E79 1,^3 FOCUS ,1211 CLEO, ESI BaR4R,E21 and BelltPI 

have set upper limits on D^-D^ mixing by analyzing the time dependence of WS 

D^ K^TT^ decays outlined in Sec. 12.2.11 Of these, the Belle limit, based on 

400 fb^^, is the most stringent. Assuming CP conservation, they find: x' < 0.72 x 

IQ-^ and -9.9 x lO'^ < y' < 6.8 x lO'^ at the 95% confidence level. 

In 2007 the BaBAR Collaboration reported evidence for D^-D^ mixing from 

(4030 ± 90) WS signal candidates and (1141500 ± 1200) RS signal candidates in 
— 1 1 2 

a 384 fb data sample, The reconstructed decay-time distribution for WS data 
and the fit results with and without mixing (assuming CP conservation) are shown 
in Fig. 1191 The fit with mixing provides a substantially better description of the 
data than the fit with no mixing. The parameters obtained from fitting the BaRAR 
data assuming CP conservation are listed in Table [3] 

In the BaBAR measurement, the significance of the mixing signal is estimated 
from the change in the log likelihood, — 2Alog£, with respect to its value at the 
global minimum. Fig. [20] shows the confidence-level contours calculated using the 
change in log likelihood from the joint estimation of two parameters. The best fit 
value of the {x' ,y') parameters to the BaBAR data is at the unphysical value of 
(x' = — 2.2x 10"'', y' = 9.7x 10"'^). As can be seen from Fig.[20l the two parameters 
are highly correlated with each other. Constraining the fit region to a:'^ > yields 
{x' = 0,y' — 6.4 X 10^"^), and corresponds to — 2Alog£ — 0.7. The no-mix point 
corresponds to — 2Alog£ = 23.9 statistical units. The maximum log likelihood is 
denoted as log£(a;' ,y'). Each systematic variation is included one at a time into 
the fit and a new log liklihood log£(a;'j , y^) is obtained. The significance of the i*'* 
systematic variation is = 2 logC{x'^ ,y') — log C{x''f , y[) /2.3, where the factor 
of 2.3 is the 68% confidence level for two degrees of freedom. Reducing — 2Alog£ 
by 1 -|- s| = 1.3 everywhere to account for systematic uncertainties results in a 
significance equivalent to 3.9 standard deviations. Predominant systematic uncer- 
tainties on the mixing parameters arise from modeling the long decay times of other 
charm decays populating the signal region and to a non-zero mean in the proper 
decay-time resolution function. 

To allow for CPV, the BaBAR analysis fits the D° and D^ samples separately to 
determine the parameters (i?J, x'^ , j/'+) and (i?^, x'^ , y'~), respectively. From 



these fitted values, the parameters i?_D = y R^Rjj and Ao = {R^—Rjj)/{R^+Rjj) 
are computed. The systematic component of the error on the BABAR measurement 
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Fig. 19. BABAR: (a) Projections of the reconstructed decay-time distribution of D" and Z?" WS 
candidates (points with error bars) and fit result integrated over the region 1.843 < rriKTr < 
1.883 GeV/c^ and 0.1445 < Am < 0.1465 GeV/c^. Assuming CP conservation, the results of fitting 
with mixing and without mixing are shown as the solid and dashed curves, respectively, (b) The 
points with error bars represent the difference between the data and the average value of the 
no-mixing fit in each data bin. The solid curve shows the difference between the fit with and 
without mixing. If there were no mixing, the data points would scatter randomly around the 
dashed horizontal line. Reprinted figure with permission from B. A uber t et al., Phys. Rev. Lett. 
98, 211802 (2007). Copyright 2007 by the American Physical Society.113 



of Ajj is mainly due to uncertainties in modeling the slight asymmetry between the 
interactions of and mesons in the detector. 

The CDF Collaboration, using a 1.5 fb~^ data sample of pp collisions at y/s = 
1.96 TeV has shown evidence for mixing in the — > K^n^ channel.El Since the 
CDF experiment was not running on the T(45) as BaBAR and Belle were, removal 
oi B ^ D decays was considerably more challenging than applying a simple center 
of mass momentum cut, as was done in the B-factory measurements. On the other 
hand, due to the much larger average boost, the average flight distance in the lab is 
greater than in the i?-factory experiments. Despite the vastly different environment, 
the central values of the mixing parameters shown in Table |3] and the (x'^, y') C.L. 
contours shown in Fig. I2H both from the CDF experiment, agree remarkably well 
with the corresponding BaBAR results. There is no evidence for CPV from any of 
the reported measurements. 
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Fig. 20. BABAR: The best fit value for the mixing parameters {x'^ , y') (•), the no-mix point (0, 0) 
(+), and confidence level (C.L.) contours evaluated for (1 — C.L.) = 0.317 (Icr), 4.55 X 10"-^ (20-), 
2.70 X 10-3 (So-), 6.33 x IQ-^ (4(7) and 5.73 x 10"'' (5o-), using the change -2Alog£ from the 
joint estimation of two parameters. The contours include the estimated systematic uncertainty. 
Reprinted figure with permission from B. Au ber t et ai, Phys. Rev. Lett. 98, 211802 (2007). Copy- 
right 2007 by the American Physical Society.HSl 
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Fig. 21. CDF: The unconstrained best fit value for the mixing parameters {x'^,y') (•), requiring 
(x'^ > 0) (0), the (0,0) (no-mix) point (-I-) and Bayesian probability contours corresponding to 
one through four equivalent Gaussian standard deviations. Reprinted figure with permission from 
T. Aaltonen et al, Phys. Rev. Lett. 100, 121802 (2008). Copyright 2008 by the American Physical 
SocietyEl 
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Table 3. D^-D" mixing and CPV parameters from D'^ — )■ JC+7r~ decays. For results with two reported 
uncertainty components, the first is statistical and the second is systematic. The results with a single 
uncertainty component include both statistical and systematic uncertainties. 



Fit type 


Parameter 


BABA^ 


Fit Results (/lO 
CDF 14 


Beiiizni 


No CPV or mixing 


Rd 


3.53 ±0.08 ±0.04 


4.15 ± 0.10 


3.77 ±0.08 ±0.05 


No CPV 


Rd 


3.03 ±0.16 ±0.10 


3.04 ±0.55 


3.64 ±0.17 


Significance 


y' 


-0.22 ±0.30 ±0.21 
9.7 ±4.4 ±3.1 
3.9 


-0.12 ±0.35 
8.5 ± 7.6 
3.8 


"■"-3.9 
2.0 


CPV allowed 


Rd 


3.03 ±0.16 ±0.10 








An 


-21 ±52 ±15 




23 ±47 




Am 






670 ± 1200 




x'++" 


-0.24 ±0.43 ±0.30 








y'+ 


9.8 ±6.4 ±4.5 








2 

x'-~ 


-0.20 ±0.41 ±0.29 








y'- 


9.6±6.1±4.3 








y' 






< 0.72 (95% C.L.) 
-28 < < 21 (95% C.L.) 



3.2.2. Lifetime Ratio Results 

E791 FOCU S, CLEO, Belle, and BaBAR have published ycp results 
| 13 | 42 | 44 | 45 | 7 2173j ^^^^ lifetime measurements ol D" K+R- and L>" ^ tt+tt- 
decays relative to that from — K~tt^ decays. Fits to the proper time distri- 
butions from the BABAR untagged analysiJ^^ are shown in Fig. [22l The lifetime 
measurement is determined from a fit using the decay time and decay-time error 
of candidates in a signal region, as described in Sec. 12.2.21 The charm background 
component shape and yield in the signal region are obtained from MC simulated 
events. The combinatorial background component shape in the signal region is esti- 
mated from sideband data. A fit to the data mass distribution is performed over 
the full mass range to estimate the total background and signal yields in the signal 
region. The combinatorial yield in the signal region is then obtained by subtracting 
the charm background yield from the total background yield there. 

The results from the Belle and BaBAR experiments on Ay show no evidence for 
CP violation. These measurements are shown at the bottom of Table ID Results for 
ycp are shown in Fig. 1231 and are giv en in T able S) The Heavy Flavor Averaging 
Group (HFAG) world average ycp 

valuePEB 

is more than four standard deviations 
away from the no- mixing hypothesis {ycp = 0). 
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Fig. 22. BABAR decay-time distributions of K' n+ (left), and K+K' (right), 

from untagged decays. In each plot, the total fit is shown as a solid line. Fit components 
are signal (white), combinatorics (gray), and charm background (black). Reprinted figures with 
permission from B. Aubert et al, Phys. Rev. D 80, 071103(R) (2009). Copyright 2009 by the 
American Physical Society."^^ 
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3.420 ±1.390 ±0.740 % 



-1.200 ±2.500 ±1.400"? 



-o.soo± 1.000 ±0.800*; 



1.310 ±0.320 ± 0.250 % 



0.110 ±0.610 ± 0.520 % 



1.160 ±0.220 ±0.180 % 



0.010 ±0.300 ±0.150 % 



0.260 ± 0.360 ±0.0e 



0.123 ± 0.248 % 



-4 -3 -2 -1 1 2 3 4 5 



-0.2 0.2 0.4 0.6 



Fig. 23. HFAG compendium of ycp (left) and A-p (right) measurements, along with world average 
values. From Refs. \lE\ [7l] 
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Table 4. R esults from tyr-P CA H and Ar measurements from E791, FOCUS, CLEO, Belle, and BABAR 
experimentsP^ ^^ I ^^ I ^^I^^ I ^*^! Measurement uncertainties are given as statistical (first) and systematic 
(second). The world-average uncertainty is statistical and systematic combined. 



Experiment 


Parameter 


Result (%) 


data sample 


E79lH2] 


AF 


0.04 ±0.14 ±0.05 


500 GeV ttN interactions (2 X 10^° events) 


FOCUS^ 


ycp 


3.42 ± 1.39 ± 0.74 


7N interactions (1 X 10'' rcconstr. D — > Kn^n) 


CLEcUiii 


ycp 


1.2 ±2.5 ±1.4 


9.0 fb~^ near T(4S') resonance ; untagged 


BelliL^ 
BellJEl 
BelliSI 


ycp 
ycp 
ycp 


0.5±1.0+;j-i 
1.31 ±0.32 ±0.25 
0.11 ±0.61 ±0.52 


23.4 fb~^ near T(45) resonance; untagged 
540 fb~^ near T(AS) resonance ; D* tagged 
673 fb-i near T{4S) resonance ; D -s- {K+K)K° 


babaM 

BABA^ 
BABA^ 


ycp 
ycp 
ycp 


1.03 ±0.33 ±0.19 
1.12 ±0.26 ±0.22 
1.16 ±0.22 ±0.18 


384 fb~^ near T{4S) resonance ; D* tagged 
384 fb~^ near Y{iS) resonance ; untagged 
D* tagged + untagged combined 


HFAd^ 


ycp 


1.107 ±0.217 


World Average 


BelleffS 


Ar 


0.01 ±0.30 ±0.15 


540 fbl near T(45) resonance ; D* tagged 


BABAI^ 


Ar 


0.26 ±0.36 ±0.08 


384 fbl near T(45) resonance ; D* tagged 


HFAd^B 


Ar 


0.123 ± 0.248 


World Average 



3.3. Results from Time- dependent Analyses of Multi-body Decays 
3.3.1. K°h+h- Analysis Results 

CLEO^I, BellfSSEH^ 

and J3aBAF0S1 have published results from time-dependent 
Dalitz-plot analyses oi D° K°TT+Tr- decays and D° K^K+K' (Belle, BaBAR). 
The CLEO Collaboration, which pioneered this technique, set 95% confidence level 
(CL) limits on the mixing parameters x and y. Results for x and y are given in 
Table [5] and world averages in Fig. [25l No evidence for CP violation has been seen 
by any of the experiments. 

Experimental systematics include variations in background and PDF models, 
efficiencies, event selection criteria, and experimental resolution effects. In addition 
to these, systematics from the chosen resonance decay model are evaluated as well. 
Additional cross-checks are performed. Fitted values of background fractions, 
lifetimes, and decay-time scale factors are determined to be consistent with ex- 
pectations or previous results. Decay-time distributions are shown in Fig. 1241 Belle, 
CLEO, and BaBAR perform additional mixing fits to check for CP violation. BaBAR 
performs separate fits to D'^ and D'^ decays, and no evidence for CP violation in 
mixing is seen. CLEO performs separate fits to and Z?" samples. Belle incorpo- 
rates additional CP-violating parameters in their mixing fit. These additional fits 
yield results consistent with the nominal fitted values. No evidence for CP violation 
is seen in any of the measurements. 
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Fig. 24. Reconstructed K^h^h decay-time distributions. Top: Belle K^tt'^tt . Bottom left: 
BABAR, ifgTT+Tr". Bottom right: BABAR, K^K^ K~ . Reprinted figures with permission from L.M. 
Zhang et al, Phys. Rev. Lett. 99, 131803 (2007) and P. del Amo Sanchez et al . Phvs . Rev. Lett. 
105, 081803 (2010) Copyright 2007 and 2010 by the American Physical Society.l^Sl^ 



3.3.2. Z?" Ktttt^ Analysis Results 

The first observation of the WS decay mode D° — > K^tt^tt^ was reported by CLEO 
in 2001.1^ Using a 9 fb~^ dataset of e+e~ collisions near the T{AS) resonance, 
they observed the decay with a 4.9 standard deviation significance and reported a 
wrong-sign rate of R^^'f = [0-43t^:iJ ± 0.07]%. Using 281 fb~^ of e+e" coUiding- 
beam data near the T{AS), Belle (2005) reported ^3 a WS branching fraction of 
Ri^Y = [0-229 ± 0.015([^o°9 ^]%. In 2006 BaBAR reported ^21 a measurement of 
R^Y" = [0-214 ± 0.008 ± 0.008]% No evidence for CP violation was observed in 
these studies. 

BaBAR analyzed the Ktttt^ decay mode using two different methods (see 
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Table 5. Results for x, y, and ycp from — > A'^tt+tt" and — > K^K~^ K~ time-dependent 
analyses. All results are from Dalitz-plot analyses except the Belle ycp result. Uncertainties 
on X and y are statistical, experimental systematic, and resonance decay model systematic, 
respectively. Prom Refs. [46] 113 HH] 



Experiment Fit Type Parameter Result 95% C.L. Limit 



CLECBZI 


No CPV 


X (%) 


l-Stl'i ±0.4 ±0.4 


(-4.7,8.6) 


CLECEZI 


No CPV 


y (%) 


-1.4^2^ ±0.8±0.4 


(-6.3,3.7) 


CLECHZI 


CPV 


X (%) 


2.3l3;4±0.4±0.4 


(-4.5,9.3) 


CLECHS 


CPV 


y (%) 


-1.5+24 ±0.8 ±0.4 


(-6.4,3.6) 


BcUe'^ 


No CPV 


X (%) 


0.80±o:29lH?toM 


(0.0, 1.6) 


Belled 


No CPV 


y (%) 


0.33±0.24lO;f2lH^ 


(-0.34,0.96) 


Belliia 


CPV 


X (%) 


0.81±0.30l°;J?in^ 


|x| < 1.6 


BellJSl 


CPV 


y (%) 


0.37±0.25tO:;^lH^ 


\y\ < 1-04 


Bell&43 


CPV 


Ii/pI 


0.86«-JgtH^±0.08 




BelliiSl 


CPV 


arg((3/p) (°) 


, ,+16+5+2 
^^-18-3-4 




BellM 


No CPV 


ycp{%) 


0.11 ±0.61 ±0.52 




BABA^ 


No CPV 


X (%) 


0.16 ±0.23 ±0.12 ±0.08 




BABAI^ 


No CPV 


y (%) 


0.57 ±0.20 ±0.13 ±0.07 




World average! '^1 ' 


No CPV 


X (%) 


0.419 ±0.211 






No CPV 


y (%) 


0.456 ± 0.186 





CLEO 2005/2007 



World average 



HFAG-cfta 


g , , , , 1 , , , , 1 , , , , 1 , , , , , , 

3 






h 




. . 1 . , . . 1 . 


. . 1 , . , . 1 . . . , 1 , . . . 1 . . . . 1 . 



1.90U + 3.300 + 0.566 % CLEO 2005/2007 



0.800 + 0.290 ± 0.170 % Bolle 2007 



0.160 + 0.230 + 0.144 % BaBar 2010 



0.419 + 0.211 % World average 



-1 1 2 3 4 5 

x(%) 




-1.400 + 2.400 + 0.894 



0.330 + 0.240 + 0.150 



0.570 + 0.200 + 0.148 



-» 0.456 + 0.186 % 



Fig. 25. K^h~^h world averages for x and y, assuming no CP violation. From Refs. [TSl 1711 



Sec. l2.3.2|) . Method I measured the time- integrated mixing rate R^i — [0.0231qq|^4± 
0.004]% with a 95% CL upper hmit of Rm < 0.054% assuming CP conservation. 
The resuh is compatible with the no-mixing hypothesis at the 4.5% CL. Additional 
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results are given in Table [6] These include measurements allowing for CP violation. 

This study estimated systematic uncertainties by varying selection cuts, chang- 
ing background PDF shapes and the decay-time resolution model, varying the D'^ 
lifetime, and changing efficiency corrections, and by performing the fit over the full 
Dalitz-plot phase space. Method I results are statistics-limited. 



Table 6. Mixing and CP violation results from BABAR analyses of 
using method I and method II (see text for description of 

parameters). 



CP conservation assumed CP violation allowed 

BABAR method ISfll 
Rm = (OmstoMi ± 0.004)% Rm = (0.010i;;;!5o? ± 0.003)% 

ay' = -0.012l^;[5g^ ± 0.002 ay' cos 4> = -0.012lg;g^f ± 0.002 

/3x'sin0 = 0.003l°:oog ± 0.000 



|p/g| = 2.2tij±0.1 



0.005 

1.0 ^ ' 



BABAR method IW^ 
0-026ltO;S ± 0.0039 x'+ = 0.0253l»;°0|^ ± 0.0039 
-0.00061°;^ ± 0.0034 y'+^^o = -O.OOOSl^Se? ± 0-0050 
^K..o = 0.0355i°;~83 ± 0.0065 
= -0.0054l0:00«± 0.0041 



■g = {5.25Io;35; ±0.12) X 10- 



Method II measured the mixing rate parameters a;'^^^o and assuming 
CP conservation, and x'^^^o and y'^T^^a allowing for CP violation, where the ) 
sign denotes measurement using only Z?'^ {D^) candidates. The method II results 
are inconsistent with the no-mixing hypothesis at a significance of 3.2 standard 
deviations. Mixing parameter values are given in Table [51 

The fit method is validated by generating Monte Carlo datasets using values for 
the PDF parameters and amplitudes taken from fits to the data. Toy Monte Carlo 
datasets are generated over the range [—0.6, 0.6] in [i^^^o/^'o, Wj^^T^o/fa], and fits to 
them reconstruct the generated values to within an offset of 30% of the statistical 
error. Studies show that increased Monte Carlo statistics reduce this offset. The 
final results include a correction for this effect. 

Systematic tests performed in the method II analysis include setting the decay- 
time resolution mean to zero (fitted value 4.2 ±0.7 fsec); changing the isobar model 
by varying the masses and widths of the included resonances within their errors; 
varying the numbers of signal and background events in each category; and changing 
the definition of the signal-region and candidate-selection criteria. 



3.3.3. D'^ — > K'^TT TT+TT Analysis Results 

BaBAR has reported a preliminary resullP^ using a time-dependent analysis of the 
D'^ K3tt mode to set a limit on the mixing rate of Rm < 0.048% at the 95% 
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confidence level (Rm — [0.019^Q Qi[g±0.002]%), and estimates that this is consistent 
with a no-mixing hypothesis at the 4.5% confidence level. This preliminary result, 
based on a 230 fb~^ data sample, used tagged events to determine the produc- 
tion flavor of the D° candidate via the charge of the slow pion tTs and reconstructed 
the candidate mass itiktvktt, the mass difference Am, the decay time t, and its 
uncertainty at- 




Fig. 26. Likelihood contours from the BABAR K3n analysis for the CP conserving fit for Ro vs. 
the mixing rate Rm = (x^ + ?/^)/2. Solid fine: A\nC = 1.15; dotted fine: AlnC = 3.0. From 
Ref.El] 

First, an unbinned extended maximum-likelihood fit is performed simultaneously 
to the RS and WS two-dimensional distributions (mi<-,rir7rj Am). Approximately 
3.5 X 10^ RS signal candidates are found for D'^ and respectively, and about 
1100 WS signal candidates for each. The much larger RS sample effectively deter- 
mines the resolution function parameters and lifetime used in the WS mixing 
fit. Backgrounds include mis-reconstructed D*^ decays, mis-reconstructed charm 
decays, and combinatorics. 

Second, fits to the decay-time distributions are performed. The wrong-sign time- 
dependence is fitted to 

- — J-- ^ Rd + ay \l RD[Tt) ^ (Ft), (71) 

where a tilde accent indicates quantities that are integrated over a region of phase 
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Fig. 27. Likelihood contours from the BABAR K3n analysis for the CP conserving fit for the 
interference term vs. the mixing rate Rm = (x^ + y'^)/2. Solid line: AlnC = 1.15; dotted line: 
A In £ = 3.0. From Ref. [75] 



space. The factor a accounts for strong-phase variation over the region; A CP- 
conserving fit (which considers and D'^ candidates simultaneously) and a fit 
that is sensitive to CP violation (which treats them separately) are performed. The 
substitutions 

cty ~> \p /ql"^ {ay cos (j)± (3 x sin (f)), (72) 
x' + y'^\p/q\^\x'+y') (73) 

are made in Eq. [Til choosing the "+" ("— ") sign for £>° (D°) candidate decays, 
respectively. The factor j3 accounts for variation over the phase space region. 

Systematic uncertainties are estimated by varying the selection criteria, the PDF 
parametrization of the decay-time resolution function, the background PDF shapes, 
and the measured Z?" lifetime value. The combined systematics are smaller than the 
statistical errors by a factor of five. 

Assuming CP conservation, the BaBAR preliminary analysis finds Rm = (x'^ + 
?/^)/2 and the effective mixing parameter ay' to be 

Rm = [0.019t°;°i^(stat.) ± 0.002(syst.)]%, (74) 
ay' = -0.006 ± 0.005(stat.) ± 0.001(syst.), (75) 



which are consistent with the no-mixing hypothesis at the 4.3% C.L. Allowing for 
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CP violation, the analysis finds 

Rm = [0.017j:°:°J^(stat.) ± 0.003(syst.)]%, (76) 

\p/q\ = l.lt^;n(stat.)±0.1(syst.), (77) 

ay'cos4> = -0.006tJ]:;]^^(stat.) ± 0.006(syst.), (78) 

cos 4> = 0.002;g ™^(stat.) ± 0.006(syst.). (79) 

Two-dimensional coverage probabilities of 68.3% and 95.0% (Alog£ = 1.15, 
3.0, respectively) are shown in Fig. [55] for the doubly Cabibbo-suppressed rate Rd 
vs. the mixing rate Rm, and in Fig. [57] for the interference term ay' / \/ +2/^ vs. 
Rm- 

3.3.4. Semileptonic Decays Analysis Results 

Results from five analyses using semileptonic decays are summarized in Ta- 
ble[7] These include results from tt-N cohisions (E79lE2at FNAL, 2 x 10^° events) 
and e+e^ interactions near the T{4S) resonance (CLEO II.V,ESl9.0 fb"^ Belle,Ell 
492 fb"^; BAa4Rsingiy-tagged,EHl87 fb"!; and BaBAR doubly-tagged,^! 344 fb"^). 
The world average is shown in Fig. [28] Results from both K'^*^ev and K'^*^ decay 
modes are included. 



Table 7. Mixing results using semileptonic D" decay modes. Uncertainties are statistical 
(first) and systematic (second), except as noted. 



Experiment 


D" modes 


Results 


E79ll62] 


Kev 
Kfiv 

Combined 


Rm = (0.16^0*?)% 
Rm = (0.0610-^^)% 
Rm = 0.llt°ofr% 
Rm < 0.50% at 90% CL 


CLEO 


Kev 

K*ev 

Combined 


Rj^l = (1.10 ±0.76)% (stat.-l-syst. combined) 
Rm = ( 0.0 ± 0.31)% (stat.-l-syst. combined) 
Rj^f = (0.16 ±0.29)% (stat.+syst. combined) 


BelldSlI 


K'-*'^ev 
K'-*^fiv 
Combined 


Rj^ = (-0.6±2.7l^-j) X 10-* 
Rm = (5.9±3.7ti^;^j x lO"" 
Rm = (1.3 ± 2.2 ± 2.0) x lO"* 
Rm < 6.1 X 10"" at 90% CL 


BaBar singly taggetJSS] 


K'-*'^ev 


Rm = 0.0023 ± 0.0012 ± 0.0004 
Rm < 0.0042 at 90% CL 


BaBar doubly taggectS^ 


K^*'>ev 


Rj^j = 0.4 X 10^* (central value) 
Rm in (-2.2,2.8) x lO"" at 68% CL 
Rm in (-13, 12) x 10-* at 90% CL 



Since the E791 analysis uses identical selections for both RS and WS candi- 
dates, systematic uncertainties in the mixing rate measurement largely cancel. Two 
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E791 1996 



CLEO 2005 



BaBar 2007 



BeUe2008 



World average 



0.110*! 



0.160 ± 0.290 ± 0.290 % 



.0.060 % 



0.013 ± 0.022 ± 0.020 % 



0.013 ± 0.027 ' 



-0.2-0.1 0.1 0.2 0.3 0.4 0.5 

Rm(%) 

Fig. 28. World average of semileptonic mixing results for Rm in Table [7] Due to possible cor- 
relations with the BABAR doubly-tagged result, the BABAR singly-tagged result is not included in 
the average. From Ref. 1151 



sources of systematic uncertainty that were investigated are the decay-time resolu- 
tion modeling and feedthrough of hadronic decays into the semileptonic sample. The 
decay-time determination is subject to detector effects and to the ambiguity from 
the missing neutrino. Decay times were estimated to be uncertain with a Gaussian 
smearing of about 15%. This affects the final mixing result by only about 10% of 
its statistical uncertainty, and is not significant. Feed-through of hadronic events 
via a hadron mis-identified as a lepton could increase the number of either RS or 
WS events. The former would overestimate the size of the RS signal and cause an 
incorrect estimate of the sensitivity to mixing, while the latter would cause a false 
WS signal. Since RS feed-through was estimated to be very small (about 3%) and, 
since no WS signal was seen, no corrections were made. Evaluation of the fit mod- 
eling systematic uncertainty (performed by adding 10 to 50 simulated mixed events 
to the WS sample) showed a bias of 10-15% toward a larger mixing rate. Since the 
final result is an upper limit, no correction was applied. 

In the Belle semileptonic result, the main systematics include uncertainty in the 
signal and background Am distributions, the amount of RS and WS backgrounds, 
RS and WS efficiencies, and modeling of the decay-time distribution. These are 
estimated separately for each of four subsamples, which are categorized by whether 
the candidate contains an electron or a muon, and which of two silicon vertex 
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detector configurations was used to record the event. The overall muon sample 
systematics are about double that of the electron samples, due in part to larger 
backgrounds in the signal regions. 

The CLEO semileptonic analysis uses simulated events to model the background 
and signal shapes used in the fit. The largest systematic comes from the statistics 
of the simulation. The second largest systematic is the shape of the decay-time dis- 
tribution, also obtained from simulation. Other sources of systematic uncertainties 
are the Q shape {Q — energy released in the D*'^ decay), electron identification, 
and fit modeling. 

The BaBAR singly-tagged semileptonic result systematics include contributions 
from signal and background PDF shapes, the decay-time resolution model, and 
decay-time PDFs for background charm decays. Other possible contributions to the 
systematic on Rm are evaluated and shown to provide no significant contribution. 
The total systematic error on Rm is about 1/3 of the statistical uncertainty. 

The BaBAR doubly-tagged semileptonic analysis finds three mixing signal can- 
didates where 2.85 background events are expected. A 50% systematic to the back- 
ground rate is assigned by comparing ten background control samples with cor- 
responding MC simulations. Other contributions to systematic uncertainties are 
ignored in comparison to this 50% error. Confidence levels are then calculated for 
Rm using a frequentist method. 

3.4. World Average Results 

From HFAG, the world average values for the D^-D^ mixing and CPV parameters 
are shown in Tab. [51 The probability contours, including both statistical and sys- 



Table 8. HFAG world average mixing and CPV parameter 
values.ED 



Parameter 


No CPV 


No direct CPV 


CPy-allowed 


X (%) 

y(%) 


U-00_o.i9 

0.74 ±0.12 


0.63 ±0.19 
0.75 ± 0.12 


Qgo+0.19 
U.DO_Q 20 

0.75 ±0.12 


\i/p\ 

fn 




1.02 ± 0.04 


0.89l«;;^ 



tematic uncertainties and allowing for CPV, are shown in Fig. [29] for the mixing 
parameters (x, y) and CPV parameters (I'z/pl, </? = aig(ci/p)). The world average {x, 
y) excludes the no-mixing point {x = Q,y ~ Q) by 10.1 standard deviations. To date, 
however, no single measurement exceeds five standard deviations. The no-CPT^ 
point (Iq'/pI = 0, (/? = 0) lies within one standard deviation of the world a verage 
{\q/p\, if) value. The recent LHCb measurements of mixing and direct (^pyEHIlZj 
are not included in these averages. 
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4. What's Next? 

Our understanding of charm physics has made great progress since 1975, when 
the first evidence for charm mesons was observed. The fact that the no-mixing 
hypothesis has been excluded by 10 standard deviations, combined with improved 
understanding of the mechanisms leading to D^-D^ mixing and CPV , leaves us in 
a position to make considerable progress in the next few years in the charm sector 
in both experimental accuracy and theoretical interpretation. Here we survey a few 
experiments that are likely to further our knowledge of charm physics and mixing 
over the next few years. 



4.1. BES-III 

With well over two years of data-taking at the time of this writing, the BES-III 
experiment at BEPC-Il has already surpassed both CLEO-c and BES-II in recorded 
luminosity. With over 200 million J/%1) and 100 million ■0' events, the experiment 
has about 4x the data samples of BES-II and CLEO-c. Performance of the machine 
is good, with peak luminosities of order 10'^^ cm~^sec~^. 

Of interest in the context of charm mixing is the machine's performance near 
the '0(3770), where it has reached a peak luminosity of 5.6 x 10^^ cm~^sec~^ and 
recorded over 1 fb~^ of data in less than a year. BES-III plans to increase the 
?/'(3770) dataset to 2.5 fb~^ in the next year or so, with a goal to eventually reach 
10 fb^^. Using the coherent decay techniques discussed earlier, it is clear that 
BEPC-II and BES-III should be able to substantially improve our knowledge of 
charm mixing in the very near future. 
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4.2. LHCb 

At the time of this writing, LHCb has embarked on its charm physics program 

using dat a take n in 2010 and 2011 and has reported results on open charm 
productioiPSEQl 

and other measurements. Given the detector design which is opti- 
mized for heavy-flavor physics, LHCb is expected to provide precision measurements 
of charm mixing and CP violation parameters in the next few years. First results 
showing evidence for direct CP violation by measuring the difference between the 
two time-integrated CP asy mmetries A{D'^ — > K^K~) and A{D^ — >■ tt+tt") have 
already been reported.^^HEl] j^esults on ycp and Ar are expected soon,fSQ] i-,agg(j on 
analyses of the higher-statistics 2011 and 2012 data samples. Additional competitive 
charm mixing and CPV measurements are expected to be forthcoming. 

4.3. SuperKEKB and Belle II 

After a decade-long successful program, the Belle detector and the KEKB acceler- 
ator stopped operations in June 2010.1^ Construction of SuperKEKB has started, 
and work has begun on the Belle II detector. An initial data sample of 5 ab~^ is 
planned to be recorded starting in 2014 with the eventual goal to reach 50 ab~^ 

lOOl 

by 2021-2022.*^ With these integrated luminosities. Belle II will have an excellent 
opportunity to improve on current D^-D^ mixing and CPV measurements. With 
5 ab~^, Belle II is expected to improve the existing statistics- limited measurements 
of X an d y by approximately a factor of two; with 50 ab~^, an additional factor of 
two.lHl 

4.4. The Super Flavor Factory SuperB 

The recently approved SuperB projectP^ in Italy will be able to contribute sub- 
stantially to our knowledge of charm mixing and CP violation. Plans call for the 
SuperB facility to be able to run at the V(3770), where a sample of 2 x 10^ D^-D^ 
pairs is expected to be accumulated. Both avenues are likely to lead to greatly in- 
creased understanding of the details of mixing and CPV . Also, like its predecessor 
JBaBAR, SuperB will be able to make use of the large charm production cross section 
near the T{AS). Estimates of statistical uncertainties using both Kit and lifetime 
ratio methods range from 6x to 12 x improvements over existing measurements, 
and possibly even better, depending on how much SuperB's improved decay-time 
resolution contributes. 

5. Summary 

Evidence for charm mixing at the level of 1% in the mixing parameters, first reported 
in 2007 by the BaBAR and Belle experiments, along with the recent evidence for 
direct CPV obtained by the LHCb Collaboration, has created renewed interest in 
the charm sector as a window to new physics. In the near future, BES-III and LHCb 
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should be reporting new charm results, along with final contributions from BaBAR, 
Belle, CDF, and CLEO. In the next several years, SupcrKEKB and SupcrB should 
improve the precision of mixing and CPV measurements by a factor of ten or more. 
This will be an exciting time for anyone interested in charm physics or precision 
flavor physics in general. 
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